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How accurate is pair additivity in describing interactions between soft polymer-based nanopar- 
ticles? Using numerical simulations we compute the free energy cost required to overlap multiple 
chains in the same region of space, and provide a quantitative measure of the effectiveness of pair 
additivity as a function of chain number and length. Our data suggest that pair additivity can 
indeed become quite inadequate as the chain density in the overlapping region increases. We also 
show that even a scaling theory based on polymer confinement can only partially account for the 
complexity of the problem. In fact, we unveil and characterize an isotropic to star-polymer cross-over 
taking place for large number of chains, and propose a revised scaling theory that better captures 
the physics of the problem. 



The question of how nano and mesoscopic particles 
spontaneously organize into complex macroscopic struc- 
tures can be considered as one of the great challenges 
in the field of soft matter today. In fact, the prospect 
of developing materials with new and exciting optical, 
mechanical, and electronic properties via the process of 
self-assembly relies on our ability to predict and control 
the phase behavior of complex fluids. 

Although most of the work on self-assembly has his- 
torically focused on small molecules, the last decade 
has witnessed several breakthroughs in particle synthe- 
sis at the meso-scale [U^ making possible the produc- 
tion of nanoparticles with complex chemical and geomet- 
rical properties with an unprecedented degree of preci- 
sion. More recently, a burst of research activities has fo- 
cused on polymer-based nanoparticles [5]. What makes 
these particles very interesting is that they break the 
dogma of mutual non-penetrability. Unlike regular col- 
loids, for which excluded volume interactions are strictly 
enforced via a hard-core or a power-law potential, com- 
plex mesoparticles such as charged or neutral star poly- 
mers, dendrimers or microgels present a very peculiar 
pair potential describing their volume interactions. Be- 
yond the small and local deformation limit, which can 
be easily described in terms of classical elasticity theory, 
these soft polymer-based nanoparticles have very unusual 
interaction potentials which allow for even complete over- 
lap among themselves [IHl]. 

Surprisingly, the simple relaxation of excluded volume 
constraints results in an extremely rich phase behavior 
that has been reported in several publications [5HT7] . No- 
tably, it was found that some classes of soft interactions 
lead to reentrant melting transitions, others to polymor- 
phic cluster phases |18j . and in general to multiple tran- 
sitions involving close-packed and non-close-packed crys- 
talline phases p^lfST] as a function of the system den- 
sity. Remarkably, the phase behavior of these systems 
is very much dependent on the the shape of the pair 
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potential, and it is today feasible to engineer interac- 
tions between star polymers or dendrimers by controlling 
their overall chemical/topological properties [22]. Likos 
et al. [23] established a criterion to predict whether for 
a bounded and repulsive potential re-entrant melting or 
cluster phases will occur based on the sign of the Fourier 
transform of the interaction. For a recent review on the 
subject we refer the reader to reference [6]. 

Given the complexity of these nanoparticles, their in- 
teractions are usually extracted via an explicit coarse- 
grained procedure to reduce the problem to a simpler 
one consisting of single effective particles interacting via 
an ad-hoc pair potential. Once a pair potential as a func- 
tion of separation r, F{2,r) is extracted, pair additivity 
among any two effective particles is assumed and phase 
diagrams are computed. 

While it is by now clear that many body effects can 
lead to density-dependent interactions between poly- 
mer pairs {25j , and methods to systematically include 
such deviations have recently been put forward ,24. ,26] 
(and references therein), in this paper we show that the 
very assumption of pair additivity can be greatly inad- 
equate in describing the interactions among polymer- 
based nanoparticles. To prove it, we compute the to- 
tal free energy cost associated with overlapping rip self- 
avoiding polymers (effective particles) of length TV and 
show that very quickly the assumed additivity of the pair 
interactions breaks down. 

If one assumes straightforward pair additivity, the free 
energy cost required to overlap rip effective particles is 
simply proportional to the number of pair interactions 

mn,) = m^ f'^'f (1) 

Alternatively, we can describe F{np) by using the scal- 
ing theory based on polymer confinement as was first sug- 
gested by Jun et. al. i27j. The main idea is that in the di- 
lute limit overlapping two chains of length N, each having 
a radius of gyration Rq ~ iV", is equivalent to confining 
a single chain of length 27V into a spherical cavity having 
a radius equal to the radius of gyration of a single chain 
i?G- The free energy cost associated with spherical con- 
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finement of a polymer of length N into a spherical cavity 
of radius R is given by f3F ~ {Rg/R)^^^^''~^^ (where 
/3 — (kBT)^-^) [28U30j . By plugging in the denominator 
R — )■ i?G and in the numerator Rq — )■ (27V)'', one obtains 
/3F(2) ~ 2^'^/(3'^~i). This equation, easily generalizable 

to Up chains as (3F{np) ^ rip'^^^^'^ ^\ states that (a) the 
free energy cost to completely overlap two chains is in- 
dependent of their length (well established result both 
theoretically and numerically [3 [5]), and (b) that over- 
lapping Hp chains has a free energy cost that does not 
grow linearly with the number of pair interactions. 

This scaling theory, that should better account for 
the density increase within the overlapping region of the 
chains, can be re-written as 



PFiup) PF{2) 1 + 



{up ~ 2) 



3!^/(3i/-l) 



(2) 



For chains in a good solvent v « 3/5, and the value of the 
exponent is close to 2.25 and is expected to cross over to 
3 for large number of polymers |57l [221 • To establish the 
correct dependence of the free energy with the number 
of chains we performed numerical simulations. 

We modeled a chain as a sequence of N spherical beads 
(monomers) of diameter a connected sequentially with 
an entropic spring of maximum extension \/2(j. The in- 
teraction between any two monomers is described via a 
hard-core potential [JT] while the entropic spring between 
consecutive monomers has the form 



/3T4(r) 



otherwise 
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The advantage of this model is that the free energy as- 
sociated with multiple overlapping chains is fully con- 
tained in their configurational entropy. To compute 
the free energy we used the thermodynamic integration 
method [3T]. The idea is to introduce a fictitious poten- 
tial 
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which acts exclusively 
different chains (z ^ j 
mass of each chain to 
dius ro = 2(7. For A = 
interacting; as A — > oo 
The free energy due to 
can then be extracted 
gral 

F[np) : 



among monomers associated to 
), and to constrain the center of 
be within a spherical shell of ra- 
the overlapping chains are not 
we recover the system of interest, 
the polymer-polymer interactions 
by performing the following inte- 
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In practice, we perform several Monte Carlo simulations 
for numerous values of A until the hard sphere limit is 
effectively reached, and compute numerically the inte- 
gral. In all our data the largest value of A was selected 




FIG. 1: Free energy cost to overlap Up chains as a 

function of the number of chains. The dashed and 
dotted lines represent respectively the straightforward 

and the scaling-based theory from Eq. [T] and Eq. [2] 
respectively. The solid line shows the revised theory in 

Eq. [6j The inset is a zoom of the data for small Up. 



to be the one for which the total energy of the system 
was decreased to a value of the order of 10~^. The small 
tail of the integral for larger values of lambda was com- 
puted by fitting the tail of the data with a power law 
and by extending the integral until the change of F{np) 
becomes clearly negligible (this typically accounts for a 
very small amount of the total free energy) . To properly 
sample V\ our simulations were run from a minimum of 
50 million to a maximum of 250 million sweeps. To im- 
prove our statistics, after every system sweep a global 
chain rotation move is also implemented. This consists 
of randomly picking a monomer and a direction in a ran- 
dom chain, and of performing a rigid rotation around a 
randomly selected axis of a small random angle of all the 
monomers connecting the selected monomer to the end 
of the chain along the selected direction. Our simula- 
tion were performed in the NVT ensemble in the dilute 
limit {V 3> 4/37ri?Q) for polymers of length iV = 64 and 
N — 128, with Up ranging from 2 to 32. 

Our results indicate that the free energy required to 
superimpose two chains equals f3F{2) = 2.20(5) for 
N = 128 and ^F(2) = 2.31(5) for N = 64. These results 
are fully consistent with previous numerical simulations 
on similar systems, and the small difference is simply due 
to finite size effects. Figure [T] presents the core results of 
our simulations and shows how Fijip) normalized by the 
pair free energy F(2) varies with the number of poly- 
mers. Lines indicating the straightforward (Eq. [T]) and 
the scaling-based (Eq. [2]) predictions are also shown as 
a reference. Our findings clearly show that Eq. [l] be- 
comes inaccurate as soon as rip > 4, Eq [2] is accurate 
up to rip < 8 polymers, yet significant and systematic 
deviations from both predictions are apparent for larger 
values of Up. Interestingly, the dependence of the free 
energy with the number of chains becomes weaker as the 
number of chains increases. 
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FIG. 2: Average radius of gyration Rq of multiple 
overlapping chains as a function of the number of chains 
n„. The line is a fit to the data. 



This result is a bit counter intuitive because one should 
expect a stronger dependence of F{np) on Up as the 
polymer concentration is increased inside the overlapping 
area. The key issue with the scaling theory is that it as- 
sumes that the size of the confining cavity (which equals 
the radius of gyration of an unconstrained polymer Rq) 
remains constant for any number of polymers. Although 
this is a good assumption for few chains, we find that 
as Up increases the average size of the confining cavity 
becomes systematically larger. This is because chains 
do swell to minimize the number of interactions. Fig- 
ure [2] shows the average size of the system, measured by 
computing the average radius of gyration of the chains 
at large values of A, (i?G("'p)|A=oo), divided by the ra- 
dius of gyration of the reference non-interacting poly- 
mers, (-Rg(«p)|a=o)- By fitting these data to a power 
law RQ{np) — a + brip, we can explicitly account for this 
correction. The revised free energy scales as 

PF^iup) ^ {n';/{a + hn'^p)f^''''''' (6) 

with a ~ 0.3, h ~ 0.7 and a ~ 0.21. As seen in Fig.[T] this 
function successfully fits most of our numerical data up 
to Up — 16. Crucially, the scaling behavior of the radius 
of gyration with the number of chains a ~ 1/5 is con- 
sistent with that of a star polymer in the swollen regime 

1/2 

{up <^ N). This suggests that a rearrangement from 
isotropically mixing chains to demixed/localized chains 
(star-polymer) may be taking place. To investigate this 
scenario, we computed the asphericity of the polymers as 
a function of Up. This is obtained by computing the in- 
ertia tensor of each polymer, and by combining the three 
eigenvalues li, I2 and into the rotational invariant pa- 
rameter [35] 

^ {h-hf + ih-hf + ih-h? . . 

2{h + h + hy ^' 

The plot of the normalized asphericity q = q{np) / q(l) as 
well as the three normalized eigenvalues U = li{np) 




FIG. 3: Average asphericity q ioT N = 64 and 
N — 128, computed according to Eq. [Tj for multiple 
overlapping chains as a function of the number of chains 
Up. The inset shows the three normalized eigenvalues of 
the shape tensor U for N — 64. 



as a function of the number of polymers Up is shown 
in Fig. [3] For a polymer in a good solvent a first order 
epsilon expansion results, in the dilute limit and for N — >■ 
00, to a value of q close to 0.415. Accurate numerical 
simulations predicted a value closer to 0.431 '33 [H]. For 
our longest chains, we find that indeed q = 0.44(2) when 
not interacting, however, as soon as polymers interact, q 
increases significantly to more than 1.6 times its original. 
This is a clear indication that the chains are not swelling 
isotropically as Up increases, but are stretching out along 
their long axis at expenses of the other two directions. 
This morphological change of the chains is suggestive of 
the fact that chain segregation may also be occurring. 

To check for chain segregation we tracked the location 
of the chains' main axes over the course of a long sim- 
ulation trajectory, and projected it over the surface of 
a sphere centered around the system. Figure |4] tracks 
down the locations of polymers' axes for Up = 2 and 
Up — 32. For the sake of clarity we show the tracks of 
only two randomly selected chains also for the case of 
Up — 32. Clearly, when Up = 2 the chain's axes can 
explore the whole spherical surface, indicating complete 
mixing of the chains. In the latter case only a small re- 
gion of the surface is explored by the two selected chains, 
clear indication of chain segregation. 

All of our data points to the fact that our system is 
effectively behaving as a star-polymer, i.e. a system of 
polymers radiating from a central core, and segregated 
into roughly conical regions. A power law fit of our data 
including only values of the free energy for Up > 15 results 

in l3F{np) oc rip^^'^"' which is indeed compatible with star 
polymers in the semi-dilute limit [35^. 

In summary, we computed the free energy cost associ- 
ated with the complete overlap of multiple self-avoiding 
chains as a function of chain number. Our data show that 
although the free energy error associated with pair addi- 
tivity of potentials between soft polymer-based nanopar- 
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FIG. 4: Chains's axial maps. The l.h.s tracks the 
location of the main axis of two chains projected onto a 

spherical surface centered around the system over the 
course of a long Monte Carlo trajectory when only two 

chains are present. The r.h.s shows the same map for 
two randomly selected chains in a system containing 32 
overlapping polymers. In the former case the polymers' 
axes perform a random walk over the spherical surface, 
indicating complete mixing of the chains, whereas in the 

latter case the chains are clearly segregated in specific 
regions. 



mer confinement of a single chain into a spherical cavity 
can better account for the free energy cost of multiple 
overlapping polymers. Finally, we have shown that de- 
mixing and chain segregation occurs when considering 
large numbers of chains, and that in this limit, the free 
energy of a star-polymer provides a better description of 
our numerical data. 

Although our study focuses on single polymer chains, 
the main message should also hold for multi-polymer- 
based particles such as star-polymers and dendrimers, 
where deviations could be much more significant even for 
small Up and could have dramatic consequences especially 
for the phase diagrams associated with systems forming 
cluster crystals where multiple particles overlap on the 
same lattice sites. It should be stressed that our results 
can be directly linked to the chain segregation process 
described by Jun et. al. [57] for polymers under confine- 
ment. In this analogy, the formation of a star-polymer 
can be understood as chain segregation with the added 
constraint of overlapping center of masses. 



tides is of the order of just a few k^T for < 4 
(making the straightforward approach in this limit, if 
not accurate, a least reasonable), such an approximation 
severely deteriorates for larger values of Up. We have also 
shown that an extended scaling theory based on poly- 
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